Range.Maps
2025-01-08
Location of files originally: https://github.com/wpetry/USTreeAtlas
library(sp)
library(raster)
library(sf)## Linking to GEOS 3.13.0, GDAL 3.10.0, PROJ 9.5.1; sf_use_s2() is TRUElibrary(maps)
library(tidyverse)## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr 1.1.4 ✔ readr 2.1.5
## ✔ forcats 1.0.0 ✔ stringr 1.5.1
## ✔ ggplot2 3.5.1 ✔ tibble 3.2.1
## ✔ lubridate 1.9.2 ✔ tidyr 1.3.0
## ✔ purrr 1.0.2## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ tidyr::extract() masks raster::extract()
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag() masks stats::lag()
## ✖ purrr::map() masks maps::map()
## ✖ dplyr::select() masks raster::select()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errorslibrary(tigris)## To enable caching of data, set `options(tigris_use_cache = TRUE)`
## in your R script or .Rprofile.options(tigris_use_cache = TRUE)Read in the GBIF data points
gbif = read.csv("../Formatted.Data/gbif.filtered.csv")Subset for Abies balsamea
ABBA.occ = gbif %>%
filter(species == "Abies balsamea") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ABBA.range = st_read("../../USTreeAtlas/shp/abiebals/")## Reading layer `abiebals' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/abiebals'
## using driver `ESRI Shapefile'
## Simple feature collection with 442 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -116.6732 ymin: 38.40429 xmax: -52.61445 ymax: 58.79789
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ABBA.range) <- 4267 maps::map(database = "state")
plot(ABBA.range$geometry, col="red", add=T)
us_states <- states(cb = TRUE)## Retrieving data for the year 2021continental_states <- us_states %>%
filter(!NAME %in% (c("Alaska","American Samoa","Guam","Commonwealth of the Northern Mariana Islands","Hawaii","United States Virgin Islands",
"Puerto Rico")))states.map = continental_states %>%
st_as_sf %>%
st_transform(st_crs(ABBA.range))ggplot(states.map) +
geom_sf()
ggplot(ABBA.range) +
geom_sf()
ggplot()+
geom_sf(data = states.map)+
geom_sf(data = ABBA.range, col = "red")
ABBA_clipped = st_intersection(ABBA.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ABBA_clipped, col = "red")+
theme_classic()
ggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ABBA_clipped, col = "red", linewidth = 2)+
geom_point(data = ABBA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ABBA.range.pdf", width = 12, height = 8)Repeat basics of steps above for all other species
Subset for Abies fraseri
ABFR.occ = gbif %>%
filter(species == "Abies fraseri") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ABFR.range = st_read("../../USTreeAtlas/shp/abiefras/")## Reading layer `abiefras' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/abiefras'
## using driver `ESRI Shapefile'
## Simple feature collection with 8 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -83.82256 ymin: 35.03492 xmax: -81.02652 ymax: 37.00807
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ABFR.range) <- 4267 ABFR_clipped = st_intersection(ABFR.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ABFR_clipped, col = "red", linewidth = 2)+
geom_point(data = ABFR.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ABFR.range.pdf", width = 12, height = 8)Subset for Acer barbatum
ACBA3.occ = gbif %>%
filter(species == "Acer barbatum") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ACBA3.range = st_read("../../USTreeAtlas/shp/acerbarb/")## Reading layer `acerbarb' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acerbarb'
## using driver `ESRI Shapefile'
## Simple feature collection with 21 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -98.45653 ymin: 28.42883 xmax: -76.26695 ymax: 37.3846
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ACBA3.range) <- 4267 ACBA3_clipped = st_intersection(ACBA3.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ACBA3_clipped, col = "red", linewidth = 2)+
geom_point(data = ACBA3.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ACBA3.range.pdf", width = 12, height = 8)Subset for Acer leucoderme
ACLE.occ = gbif %>%
filter(species == "Acer leucoderme") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ACLE.range = st_read("../../USTreeAtlas/shp/acerleuc/")## Reading layer `acerleuc' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acerleuc'
## using driver `ESRI Shapefile'
## Simple feature collection with 21 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -95.10542 ymin: 30.2384 xmax: -78.57996 ymax: 36.43029
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ACLE.range) <- 4267 ACLE_clipped = st_intersection(ACLE.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ACLE_clipped, col = "red", linewidth = 2)+
geom_point(data = ACLE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ACLE.range.pdf", width = 12, height = 8)Subset for Acer negundo
ACNE3.occ = gbif %>%
filter(species == "Acer negundo") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ACNE3.range = st_read("../../USTreeAtlas/shp/acernegu/")## Reading layer `acernegu' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acernegu'
## using driver `ESRI Shapefile'
## Simple feature collection with 114 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -123.3583 ymin: 28.30204 xmax: -67.68153 ymax: 54.28115
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ACNE3.range) <- 4267 ACNE3_clipped = st_intersection(ACNE3.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ACNE3_clipped, col = "red", linewidth = 2)+
geom_point(data = ACNE3.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ACNE3.range.pdf", width = 12, height = 8)Subset for Acer nigrum
ACNI5.occ = gbif %>%
filter(species == "Acer nigrum") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ACNI5.range = st_read("../../USTreeAtlas/shp/acernigr/")## Reading layer `acernigr' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acernigr'
## using driver `ESRI Shapefile'
## Simple feature collection with 92 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.94963 ymin: 35.15816 xmax: -72.13584 ymax: 45.95228
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ACNI5.range) <- 4267 ACNI5_clipped = st_intersection(ACNI5.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ACNI5_clipped, col = "red", linewidth = 2)+
geom_point(data = ACNI5.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ACNI5.range.pdf", width = 12, height = 8)Subset for Acer pensylvanicum
ACPE.occ = gbif %>%
filter(species == "Acer pensylvanicum") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ACPE.range = st_read("../../USTreeAtlas/shp/acerpens/")## Reading layer `acerpens' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acerpens'
## using driver `ESRI Shapefile'
## Simple feature collection with 139 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -89.90945 ymin: 34.813 xmax: -59.68972 ymax: 49.26055
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ACPE.range) <- 4267 ACPE_clipped = st_intersection(ACPE.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ACPE_clipped, col = "red", linewidth = 2)+
geom_point(data = ACPE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ACPE.range.pdf", width = 12, height = 8)Subset for Acer rubrum
ACRU.occ = gbif %>%
filter(species == "Acer rubrum") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ACRU.range = st_read("../../USTreeAtlas/shp/acerrubr/")## Reading layer `acerrubr' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acerrubr'
## using driver `ESRI Shapefile'
## Simple feature collection with 498 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -97.25412 ymin: 25.5467 xmax: -52.61445 ymax: 49.401
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ACRU.range) <- 4267 ACRU_clipped = st_intersection(ACRU.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ACRU_clipped, col = "red", linewidth = 2)+
geom_point(data = ACRU.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ACRU.range.pdf", width = 12, height = 8)Subset for Acer saccharinum
ACSA2.occ = gbif %>%
filter(species == "Acer saccharinum") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ACSA2.range = st_read("../../USTreeAtlas/shp/acersacc/")## Reading layer `acersacc' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acersacc'
## using driver `ESRI Shapefile'
## Simple feature collection with 38 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -98.05198 ymin: 29.9023 xmax: -66.47555 ymax: 48.20345
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ACSA2.range) <- 4267 ACSA2_clipped = st_intersection(ACSA2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ACSA2_clipped, col = "red", linewidth = 2)+
geom_point(data = ACSA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ACSA2.range.pdf", width = 12, height = 8)Subset for Acer saccharum
ACSA3.occ = gbif %>%
filter(species == "Acer saccharum") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ACSA3.range = st_read("../../USTreeAtlas/shp/acersacr/")## Reading layer `acersacr' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acersacr'
## using driver `ESRI Shapefile'
## Simple feature collection with 441 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -97.7326 ymin: 33.45683 xmax: -59.68972 ymax: 49.4221
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ACSA3.range) <- 4267 ACSA3_clipped = st_intersection(ACSA3.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ACSA3_clipped, col = "red", linewidth = 2)+
geom_point(data = ACSA3.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ACSA3.range.pdf", width = 12, height = 8)Subset for Acer spicatum
ACSP2.occ = gbif %>%
filter(species == "Acer spicatum") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ACSP2.range = st_read("../../USTreeAtlas/shp/acerspic/")## Reading layer `acerspic' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acerspic'
## using driver `ESRI Shapefile'
## Simple feature collection with 264 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -104.7974 ymin: 34.86549 xmax: -52.61445 ymax: 54.83275
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ACSP2.range) <- 4267 ACSP2_clipped = st_intersection(ACSP2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ACSP2_clipped, col = "red", linewidth = 2)+
geom_point(data = ACSP2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ACSP2.range.pdf", width = 12, height = 8)Subset for Aesculus flava. Range map not available.
AEFL.occ = gbif %>%
filter(species == "Aesculus flava") %>%
select(species,decimalLatitude,decimalLongitude)ggplot()+
geom_sf(data = states.map, fill = "white")+
geom_point(data = AEFL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/AEFL.range.pdf", width = 12, height = 8)Subset for Aesculus glabra
AEGL.occ = gbif %>%
filter(species == "Aesculus glabra") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
AEGL.range = st_read("../../USTreeAtlas/shp/aescglab/")## Reading layer `aescglab' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/aescglab'
## using driver `ESRI Shapefile'
## Simple feature collection with 17 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -99.97602 ymin: 29.33802 xmax: -79.56816 ymax: 42.69205
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(AEGL.range) <- 4267 AEGL_clipped = st_intersection(AEGL.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = AEGL_clipped, col = "red", linewidth = 2)+
geom_point(data = AEGL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/AEGL.range.pdf", width = 12, height = 8)Subset for Ailanthus altissima. Range map not available.
AIAL.occ = gbif %>%
filter(species == "Ailanthus altissima") %>%
select(species,decimalLatitude,decimalLongitude)ggplot()+
geom_sf(data = states.map, fill = "white")+
geom_point(data = AIAL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/AIAL.range.pdf", width = 12, height = 8)Subset for Annona glabra. Range map not available.
ANGL4.occ = gbif %>%
filter(species == "Annona glabra") %>%
select(species,decimalLatitude,decimalLongitude)ggplot()+
geom_sf(data = states.map, fill = "white")+
geom_point(data = ANGL4.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ANGL4.range.pdf", width = 12, height = 8)Subset for Asimina triloba
ASTR.occ = gbif %>%
filter(species == "Asimina triloba") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ASTR.range = st_read("../../USTreeAtlas/shp/asimtril/")## Reading layer `asimtril' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/asimtril'
## using driver `ESRI Shapefile'
## Simple feature collection with 44 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -97.22963 ymin: 29.65825 xmax: -74.09969 ymax: 43.43197
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ASTR.range) <- 4267 ASTR_clipped = st_intersection(ASTR.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ASTR_clipped, col = "red", linewidth = 2)+
geom_point(data = ASTR.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ASTR.range.pdf", width = 12, height = 8)Subset for Betula alleghaniensis
BEAL2.occ = gbif %>%
filter(species == "Betula alleghaniensis") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
BEAL2.range = st_read("../../USTreeAtlas/shp/betualle/")## Reading layer `betualle' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/betualle'
## using driver `ESRI Shapefile'
## Simple feature collection with 91 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.07653 ymin: 34.77789 xmax: -52.61445 ymax: 49.95027
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(BEAL2.range) <- 4267 BEAL2_clipped = st_intersection(BEAL2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = BEAL2_clipped, col = "red", linewidth = 2)+
geom_point(data = BEAL2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/BEAL2.range.pdf", width = 12, height = 8)Subset for Betula lenta
BELE.occ = gbif %>%
filter(species == "Betula lenta") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
BELE.range = st_read("../../USTreeAtlas/shp/betulent/")## Reading layer `betulent' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/betulent'
## using driver `ESRI Shapefile'
## Simple feature collection with 20 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -87.99353 ymin: 33.1902 xmax: -68.21455 ymax: 48.70293
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(BELE.range) <- 4267 BELE_clipped = st_intersection(BELE.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = BELE_clipped, col = "red", linewidth = 2)+
geom_point(data = BELE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/BELE.range.pdf", width = 12, height = 8)Subset for Betula nigra
BENI.occ = gbif %>%
filter(species == "Betula nigra") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
BENI.range = st_read("../../USTreeAtlas/shp/betunigr/")## Reading layer `betunigr' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/betunigr'
## using driver `ESRI Shapefile'
## Simple feature collection with 48 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.42217 ymin: 29.56167 xmax: -70.57861 ymax: 45.09521
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(BENI.range) <- 4267 BENI_clipped = st_intersection(BENI.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = BENI_clipped, col = "red", linewidth = 2)+
geom_point(data = BENI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/BENI.range.pdf", width = 12, height = 8)Subset for Betula papyrifera
BEPA.occ = gbif %>%
filter(species == "Betula papyrifera") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
BEPA.range = st_read("../../USTreeAtlas/shp/betupapy/")## Reading layer `betupapy' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/betupapy'
## using driver `ESRI Shapefile'
## Simple feature collection with 482 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -164.7386 ymin: 35.027 xmax: -52.61445 ymax: 68.50498
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(BEPA.range) <- 4267 BEPA_clipped = st_intersection(BEPA.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = BEPA_clipped, col = "red", linewidth = 2)+
geom_point(data = BEPA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/BEPA.range.pdf", width = 12, height = 8)Subset for Betula populifolia
BEPO.occ = gbif %>%
filter(species == "Betula populifolia") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
BEPO.range = st_read("../../USTreeAtlas/shp/betupopu/")## Reading layer `betupopu' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/betupopu'
## using driver `ESRI Shapefile'
## Simple feature collection with 145 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -87.36114 ymin: 35.00822 xmax: -59.68972 ymax: 47.98228
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(BEPO.range) <- 4267 BEPO_clipped = st_intersection(BEPO.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = BEPO_clipped, col = "red", linewidth = 2)+
geom_point(data = BEPO.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/BEPO.range.pdf", width = 12, height = 8)Subset for Carpinus caroliniana
CACA18.occ = gbif %>%
filter(species == "Carpinus caroliniana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
CACA18.range = st_read("../../USTreeAtlas/shp/carpcaro/")## Reading layer `carpcaro' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/carpcaro'
## using driver `ESRI Shapefile'
## Simple feature collection with 74 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -108.0685 ymin: 13.74625 xmax: -68.47728 ymax: 47.47089
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(CACA18.range) <- 4267 CACA18_clipped = st_intersection(CACA18.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = CACA18_clipped, col = "red", linewidth = 2)+
geom_point(data = CACA18.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/CACA18.range.pdf", width = 12, height = 8)Subset for Carya alba. No range map available
CAAL27.occ = gbif %>%
filter(species == "Carya alba") %>%
select(species,decimalLatitude,decimalLongitude)ggplot()+
geom_sf(data = states.map, fill = "white")+
geom_point(data = CAAL27.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/CAAL27.range.pdf", width = 12, height = 8)Subset for Carya aquatica
CAAQ2.occ = gbif %>%
filter(species == "Carya aquatica") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
CAAQ2.range = st_read("../../USTreeAtlas/shp/caryaqua/")## Reading layer `caryaqua' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/caryaqua'
## using driver `ESRI Shapefile'
## Simple feature collection with 46 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -97.70158 ymin: 26.41639 xmax: -75.71973 ymax: 37.88053
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(CAAQ2.range) <- 4267 CAAQ2_clipped = st_intersection(CAAQ2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = CAAQ2_clipped, col = "red", linewidth = 2)+
geom_point(data = CAAQ2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/CAAQ2.range.pdf", width = 12, height = 8)Subset for Carya cordiformis
CACO15.occ = gbif %>%
filter(species == "Carya cordiformis") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
CACO15.range = st_read("../../USTreeAtlas/shp/carycord/")## Reading layer `carycord' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/carycord'
## using driver `ESRI Shapefile'
## Simple feature collection with 62 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -98.77852 ymin: 29.61891 xmax: -70.9063 ymax: 47.43512
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(CACO15.range) <- 4267 CACO15_clipped = st_intersection(CACO15.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = CACO15_clipped, col = "red", linewidth = 2)+
geom_point(data = CACO15.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/CACO15.range.pdf", width = 12, height = 8)Subset for Carya glabra
CAGL8.occ = gbif %>%
filter(species == "Carya glabra") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
CAGL8.range = st_read("../../USTreeAtlas/shp/caryglab/")## Reading layer `caryglab' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/caryglab'
## using driver `ESRI Shapefile'
## Simple feature collection with 88 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -95.0925 ymin: 26.5667 xmax: -69.92722 ymax: 44.09084
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(CAGL8.range) <- 4267 CAGL8_clipped = st_intersection(CAGL8.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = CAGL8_clipped, col = "red", linewidth = 2)+
geom_point(data = CAGL8.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/CAGL8.range.pdf", width = 12, height = 8)Subset for Carya illinoinensis
CAIL2.occ = gbif %>%
filter(species == "Carya illinoinensis") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
CAIL2.range = st_read("../../USTreeAtlas/shp/caryilli/")## Reading layer `caryilli' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/caryilli'
## using driver `ESRI Shapefile'
## Simple feature collection with 22 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -104.0874 ymin: 25.38815 xmax: -84.48216 ymax: 42.30187
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(CAIL2.range) <- 4267 CAIL2_clipped = st_intersection(CAIL2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = CAIL2_clipped, col = "red", linewidth = 2)+
geom_point(data = CAIL2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/CAIL2.range.pdf", width = 12, height = 8)Subset for Carya ovata
CAOV2.occ = gbif %>%
filter(species == "Carya ovata") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
CAOV2.range = st_read("../../USTreeAtlas/shp/caryovat/")## Reading layer `caryovat' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/caryovat'
## using driver `ESRI Shapefile'
## Simple feature collection with 49 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -101.1183 ymin: 19.93666 xmax: -69.77362 ymax: 46.39629
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(CAOV2.range) <- 4267 CAOV2_clipped = st_intersection(CAOV2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = CAOV2_clipped, col = "red", linewidth = 2)+
geom_point(data = CAOV2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/CAOV2.range.pdf", width = 12, height = 8)Subset for Carya texana
CATE9.occ = gbif %>%
filter(species == "Carya texana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
CATE9.range = st_read("../../USTreeAtlas/shp/carytexa/")## Reading layer `carytexa' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/carytexa'
## using driver `ESRI Shapefile'
## Simple feature collection with 10 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -99.2131 ymin: 28.6377 xmax: -87.35191 ymax: 40.45395
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(CATE9.range) <- 4267 CATE9_clipped = st_intersection(CATE9.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = CATE9_clipped, col = "red", linewidth = 2)+
geom_point(data = CATE9.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/CATE9.range.pdf", width = 12, height = 8)Subset for Castanea dentata
CADE12.occ = gbif %>%
filter(species == "Castanea dentata") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
CADE12.range = st_read("../../USTreeAtlas/shp/castdent/")## Reading layer `castdent' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/castdent'
## using driver `ESRI Shapefile'
## Simple feature collection with 69 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -90.06321 ymin: 30.78233 xmax: -68.80517 ymax: 44.94212
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(CADE12.range) <- 4267 CADE12_clipped = st_intersection(CADE12.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = CADE12_clipped, col = "red", linewidth = 2)+
geom_point(data = CADE12.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/CADE12.range.pdf", width = 12, height = 8)Subset for Celtis laevigata
CELA.occ = gbif %>%
filter(species == "Celtis laevigata") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
CELA.range = st_read("../../USTreeAtlas/shp/celtlaev/")## Reading layer `celtlaev' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/celtlaev'
## using driver `ESRI Shapefile'
## Simple feature collection with 83 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -104.1844 ymin: 21.94096 xmax: -75.44749 ymax: 40.79288
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(CELA.range) <- 4267 CELA_clipped = st_intersection(CELA.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = CELA_clipped, col = "red", linewidth = 2)+
geom_point(data = CELA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/CELA.range.pdf", width = 12, height = 8)Subset for Celtis occidentalis
CEOC.occ = gbif %>%
filter(species == "Celtis occidentalis") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
CEOC.range = st_read("../../USTreeAtlas/shp/celtocci/")## Reading layer `celtocci' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/celtocci'
## using driver `ESRI Shapefile'
## Simple feature collection with 63 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -104.591 ymin: 32.21578 xmax: -71.17735 ymax: 50.1872
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(CEOC.range) <- 4267 CEOC_clipped = st_intersection(CEOC.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = CEOC_clipped, col = "red", linewidth = 2)+
geom_point(data = CEOC.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/CEOC.range.pdf", width = 12, height = 8)Subset for Cercis canadensis. Range map unavailable
CECA4.occ = gbif %>%
filter(species == "Cercis canadensis") %>%
select(species,decimalLatitude,decimalLongitude)ggplot()+
geom_sf(data = states.map, fill = "white")+
geom_point(data = CECA4.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/CECA4.range.pdf", width = 12, height = 8)Subset for Chamaecyparis thyoides
CHTH2.occ = gbif %>%
filter(species == "Chamaecyparis thyoides") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
CHTH2.range = st_read("../../USTreeAtlas/shp/chamthyo/")## Reading layer `chamthyo' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/chamthyo'
## using driver `ESRI Shapefile'
## Simple feature collection with 67 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -89.87479 ymin: 29.12146 xmax: -68.76874 ymax: 44.42347
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(CHTH2.range) <- 4267 CHTH2_clipped = st_intersection(CHTH2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = CHTH2_clipped, col = "red", linewidth = 2)+
geom_point(data = CHTH2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/CHTH2.range.pdf", width = 12, height = 8)Subset for Cornus florida
COFL2.occ = gbif %>%
filter(species == "Cornus florida") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
COFL2.range = st_read("../../USTreeAtlas/shp/cornflor/")## Reading layer `cornflor' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/cornflor'
## using driver `ESRI Shapefile'
## Simple feature collection with 84 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -100.8646 ymin: 18.52992 xmax: -69.92722 ymax: 43.94442
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(COFL2.range) <- 4267 COFL2_clipped = st_intersection(COFL2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = COFL2_clipped, col = "red", linewidth = 2)+
geom_point(data = COFL2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/COFL2.range.pdf", width = 12, height = 8)Subset for Diospyros virginiana
DIVI5.occ = gbif %>%
filter(species == "Diospyros virginiana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
DIVI5.range = st_read("../../USTreeAtlas/shp/diosvirg/")## Reading layer `diosvirg' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/diosvirg'
## using driver `ESRI Shapefile'
## Simple feature collection with 100 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -99.25706 ymin: 25.02361 xmax: -72.77579 ymax: 41.37863
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(DIVI5.range) <- 4267 DIVI5_clipped = st_intersection(DIVI5.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = DIVI5_clipped, col = "red", linewidth = 2)+
geom_point(data = DIVI5.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/DIVI5.range.pdf", width = 12, height = 8)Subset for Fagus grandifolia
FAGR.occ = gbif %>%
filter(species == "Fagus grandifolia") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
FAGR.range = st_read("../../USTreeAtlas/shp/fagugran/")## Reading layer `fagugran' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/fagugran'
## using driver `ESRI Shapefile'
## Simple feature collection with 171 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -99.64175 ymin: 19.44994 xmax: -59.68972 ymax: 47.92627
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(FAGR.range) <- 4267 FAGR_clipped = st_intersection(FAGR.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = FAGR_clipped, col = "red", linewidth = 2)+
geom_point(data = FAGR.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/FAGR.range.pdf", width = 12, height = 8)Subset for Fraxinus americana
FRAM2.occ = gbif %>%
filter(species == "Fraxinus americana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
FRAM2.range = st_read("../../USTreeAtlas/shp/fraxamer/")## Reading layer `fraxamer' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/fraxamer'
## using driver `ESRI Shapefile'
## Simple feature collection with 157 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -98.82271 ymin: 28.62632 xmax: -59.79973 ymax: 48.1909
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(FRAM2.range) <- 4267 FRAM2_clipped = st_intersection(FRAM2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = FRAM2_clipped, col = "red", linewidth = 2)+
geom_point(data = FRAM2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/FRAM2.range.pdf", width = 12, height = 8)Subset for Fraxinus caroliniana
FRCA3.occ = gbif %>%
filter(species == "Fraxinus caroliniana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
FRCA3.range = st_read("../../USTreeAtlas/shp/fraxcaro/")## Reading layer `fraxcaro' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/fraxcaro'
## using driver `ESRI Shapefile'
## Simple feature collection with 71 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -95.62966 ymin: 25.74786 xmax: -75.44749 ymax: 38.66853
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(FRCA3.range) <- 4267 FRCA3_clipped = st_intersection(FRCA3.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = FRCA3_clipped, col = "red", linewidth = 2)+
geom_point(data = FRCA3.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/FRCA3.range.pdf", width = 12, height = 8)Subset for Fraxinus nigra
FRNI.occ = gbif %>%
filter(species == "Fraxinus nigra") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
FRNI.range = st_read("../../USTreeAtlas/shp/fraxnigr/")## Reading layer `fraxnigr' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/fraxnigr'
## using driver `ESRI Shapefile'
## Simple feature collection with 144 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -97.86788 ymin: 37.5547 xmax: -57.8979 ymax: 51.50689
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(FRNI.range) <- 4267 FRNI_clipped = st_intersection(FRNI.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = FRNI_clipped, col = "red", linewidth = 2)+
geom_point(data = FRNI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/FRNI.range.pdf", width = 12, height = 8)Subset for Fraxinus pennsylvanica
FRPE.occ = gbif %>%
filter(species == "Fraxinus pennsylvanica") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
FRPE.range = st_read("../../USTreeAtlas/shp/fraxpenn/")## Reading layer `fraxpenn' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/fraxpenn'
## using driver `ESRI Shapefile'
## Simple feature collection with 188 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -111.6724 ymin: 28.31472 xmax: -59.68972 ymax: 54.27539
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(FRPE.range) <- 4267 FRPE_clipped = st_intersection(FRPE.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = FRPE_clipped, col = "red", linewidth = 2)+
geom_point(data = FRPE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/FRPE.range.pdf", width = 12, height = 8)Subset for Fraxinus profunda
FRPR.occ = gbif %>%
filter(species == "Fraxinus profunda") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
FRPR.range = st_read("../../USTreeAtlas/shp/fraxprof/")## Reading layer `fraxprof' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/fraxprof'
## using driver `ESRI Shapefile'
## Simple feature collection with 58 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -93.60954 ymin: 28.59934 xmax: -75.44749 ymax: 41.38225
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(FRPR.range) <- 4267 FRPR_clipped = st_intersection(FRPR.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = FRPR_clipped, col = "red", linewidth = 2)+
geom_point(data = FRPR.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/FRPR.range.pdf", width = 12, height = 8)Subset for Fraxinus quadrangulata
FRQU.occ = gbif %>%
filter(species == "Fraxinus quadrangulata") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
FRQU.range = st_read("../../USTreeAtlas/shp/fraxquad/")## Reading layer `fraxquad' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/fraxquad'
## using driver `ESRI Shapefile'
## Simple feature collection with 18 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.2248 ymin: 33.2491 xmax: -81.49872 ymax: 42.96199
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(FRQU.range) <- 4267 FRQU_clipped = st_intersection(FRQU.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = FRQU_clipped, col = "red", linewidth = 2)+
geom_point(data = FRQU.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/FRQU.range.pdf", width = 12, height = 8)Subset for Fraxinus texensis/Fraxinus albicans (gbif)
FRTE.occ = gbif %>%
filter(species == "Fraxinus albicans") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
FRTE.range = st_read("../../USTreeAtlas/shp/fraxtexe/")## Reading layer `fraxtexe' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/fraxtexe'
## using driver `ESRI Shapefile'
## Simple feature collection with 8 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -101.2223 ymin: 28.28765 xmax: -96.66921 ymax: 34.50977
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(FRTE.range) <- 4267 FRTE_clipped = st_intersection(FRTE.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = FRTE_clipped, col = "red", linewidth = 2)+
geom_point(data = FRTE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/FRTE.range.pdf", width = 12, height = 8)Subset for Gleditsia aquatica
GLAQ.occ = gbif %>%
filter(species == "Gleditsia aquatica") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
GLAQ.range = st_read("../../USTreeAtlas/shp/gledaqua/")## Reading layer `gledaqua' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/gledaqua'
## using driver `ESRI Shapefile'
## Simple feature collection with 20 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.36879 ymin: 27.27744 xmax: -79.58623 ymax: 38.97958
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(GLAQ.range) <- 4267 GLAQ_clipped = st_intersection(GLAQ.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = GLAQ_clipped, col = "red", linewidth = 2)+
geom_point(data = GLAQ.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/GLAQ.range.pdf", width = 12, height = 8)Subset for Gleditsia triacanthos
GLTR.occ = gbif %>%
filter(species == "Gleditsia triacanthos") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
GLTR.range = st_read("../../USTreeAtlas/shp/gledtria/")## Reading layer `gledtria' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/gledtria'
## using driver `ESRI Shapefile'
## Simple feature collection with 18 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -101.5966 ymin: 28.73964 xmax: -76.56055 ymax: 43.67306
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(GLTR.range) <- 4267 GLTR_clipped = st_intersection(GLTR.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = GLTR_clipped, col = "red", linewidth = 2)+
geom_point(data = GLTR.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/GLTR.range.pdf", width = 12, height = 8)Subset for Gordonia lasianthus
GOLA.occ = gbif %>%
filter(species == "Gordonia lasianthus") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
GOLA.range = st_read("../../USTreeAtlas/shp/gordlasi/")## Reading layer `gordlasi' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/gordlasi'
## using driver `ESRI Shapefile'
## Simple feature collection with 29 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -89.66769 ymin: 27.00125 xmax: -75.71973 ymax: 35.99722
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(GOLA.range) <- 4267 GOLA_clipped = st_intersection(GOLA.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = GOLA_clipped, col = "red", linewidth = 2)+
geom_point(data = GOLA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/GOLA.range.pdf", width = 12, height = 8)Subset for Gymnocladus dioicus
GYDI.occ = gbif %>%
filter(species == "Gymnocladus dioicus") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
GYDI.range = st_read("../../USTreeAtlas/shp/gymndioi/")## Reading layer `gymndioi' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/gymndioi'
## using driver `ESRI Shapefile'
## Simple feature collection with 44 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -99.97971 ymin: 33.98936 xmax: -76.18254 ymax: 45.07224
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(GYDI.range) <- 4267 GYDI_clipped = st_intersection(GYDI.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = GYDI_clipped, col = "red", linewidth = 2)+
geom_point(data = GYDI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/GYDI.range.pdf", width = 12, height = 8)Subset for Ilex opaca
ILOP.occ = gbif %>%
filter(species == "Ilex opaca") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ILOP.range = st_read("../../USTreeAtlas/shp/ilexopac/")## Reading layer `ilexopac' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/ilexopac'
## using driver `ESRI Shapefile'
## Simple feature collection with 67 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -98.05645 ymin: 27.40893 xmax: -70.25293 ymax: 42.31931
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ILOP.range) <- 4267 ILOP_clipped = st_intersection(ILOP.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ILOP_clipped, col = "red", linewidth = 2)+
geom_point(data = ILOP.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ILOP.range.pdf", width = 12, height = 8)Subset for Juglans cinerea
JUCI.occ = gbif %>%
filter(species == "Juglans cinerea") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
JUCI.range = st_read("../../USTreeAtlas/shp/juglcine/")## Reading layer `juglcine' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/juglcine'
## using driver `ESRI Shapefile'
## Simple feature collection with 44 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -94.84976 ymin: 33.36165 xmax: -65.67808 ymax: 47.24971
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(JUCI.range) <- 4267 JUCI_clipped = st_intersection(JUCI.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = JUCI_clipped, col = "red", linewidth = 2)+
geom_point(data = JUCI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/JUCI.range.pdf", width = 12, height = 8)Subset for Juglans nigra
JUNI.occ = gbif %>%
filter(species == "Juglans nigra") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
JUNI.range = st_read("../../USTreeAtlas/shp/juglnigr/")## Reading layer `juglnigr' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/juglnigr'
## using driver `ESRI Shapefile'
## Simple feature collection with 61 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -100.5016 ymin: 28.65753 xmax: -71.5475 ymax: 45.06173
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(JUNI.range) <- 4267 JUNI_clipped = st_intersection(JUNI.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = JUNI_clipped, col = "red", linewidth = 2)+
geom_point(data = JUNI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/JUNI.range.pdf", width = 12, height = 8)Subset for Juniperus virginiana
JUVI.occ = gbif %>%
filter(species == "Juniperus virginiana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
JUVI.range = st_read("../../USTreeAtlas/shp/junivirg/")## Reading layer `junivirg' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/junivirg'
## using driver `ESRI Shapefile'
## Simple feature collection with 60 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -103.5662 ymin: 29.45831 xmax: -69.77362 ymax: 47.07206
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(JUVI.range) <- 4267 JUVI_clipped = st_intersection(JUVI.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = JUVI_clipped, col = "red", linewidth = 2)+
geom_point(data = JUVI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/JUVI.range.pdf", width = 12, height = 8)Subset for Larix laricina
LALA.occ = gbif %>%
filter(species == "Larix laricina") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
LALA.range = st_read("../../USTreeAtlas/shp/larilari/")## Reading layer `larilari' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/larilari'
## using driver `ESRI Shapefile'
## Simple feature collection with 471 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -160.2373 ymin: 39.14292 xmax: -52.61445 ymax: 67.67176
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(LALA.range) <- 4267 LALA_clipped = st_intersection(LALA.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = LALA_clipped, col = "red", linewidth = 2)+
geom_point(data = LALA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/LALA.range.pdf", width = 12, height = 8)Subset for Liquidambar styraciflua
LIST2.occ = gbif %>%
filter(species == "Liquidambar styraciflua") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
LIST2.range = st_read("../../USTreeAtlas/shp/liqustyr/")## Reading layer `liqustyr' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/liqustyr'
## using driver `ESRI Shapefile'
## Simple feature collection with 81 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -99.09172 ymin: 13.20806 xmax: -73.26896 ymax: 41.15219
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(LIST2.range) <- 4267 LIST2_clipped = st_intersection(LIST2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = LIST2_clipped, col = "red", linewidth = 2)+
geom_point(data = LIST2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/LIST2.range.pdf", width = 12, height = 8)Subset for Liriodendron tulipifera
LITU.occ = gbif %>%
filter(species == "Liriodendron tulipifera") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
LITU.range = st_read("../../USTreeAtlas/shp/lirituli/")## Reading layer `lirituli' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/lirituli'
## using driver `ESRI Shapefile'
## Simple feature collection with 61 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -93.14224 ymin: 28.39825 xmax: -71.36444 ymax: 44.53041
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(LITU.range) <- 4267 LITU_clipped = st_intersection(LITU.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = LITU_clipped, col = "red", linewidth = 2)+
geom_point(data = LITU.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/LITU.range.pdf", width = 12, height = 8)Subset for Maclura pomifera
MAPO.occ = gbif %>%
filter(species == "Maclura pomifera") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
MAPO.range = st_read("../../USTreeAtlas/shp/maclpomi/")## Reading layer `maclpomi' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/maclpomi'
## using driver `ESRI Shapefile'
## Simple feature collection with 2 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -103.3863 ymin: 29.07236 xmax: -93.57389 ymax: 34.64141
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(MAPO.range) <- 4267 MAPO_clipped = st_intersection(MAPO.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = MAPO_clipped, col = "red", linewidth = 2)+
geom_point(data = MAPO.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/MAPO.range.pdf", width = 12, height = 8)Subset for Magnolia acuminata
MAAC.occ = gbif %>%
filter(species == "Magnolia acuminata") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
MAAC.range = st_read("../../USTreeAtlas/shp/magnacum/")## Reading layer `magnacum' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/magnacum'
## using driver `ESRI Shapefile'
## Simple feature collection with 49 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -94.6578 ymin: 30.68451 xmax: -75.45156 ymax: 43.38474
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(MAAC.range) <- 4267 MAAC_clipped = st_intersection(MAAC.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = MAAC_clipped, col = "red", linewidth = 2)+
geom_point(data = MAAC.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/MAAC.range.pdf", width = 12, height = 8)Subset for Magnolia fraseri
MAFR.occ = gbif %>%
filter(species == "Magnolia fraseri") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
MAFR.range = st_read("../../USTreeAtlas/shp/magnfras/")## Reading layer `magnfras' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/magnfras'
## using driver `ESRI Shapefile'
## Simple feature collection with 4 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -84.82449 ymin: 34.46459 xmax: -77.81695 ymax: 39.62026
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(MAFR.range) <- 4267 MAFR_clipped = st_intersection(MAFR.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = MAFR_clipped, col = "red", linewidth = 2)+
geom_point(data = MAFR.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/MAFR.range.pdf", width = 12, height = 8)Subset for Magnolia grandiflora
MAGR4.occ = gbif %>%
filter(species == "Magnolia grandiflora") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
MAGR4.range = st_read("../../USTreeAtlas/shp/magngran/")## Reading layer `magngran' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/magngran'
## using driver `ESRI Shapefile'
## Simple feature collection with 28 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.0386 ymin: 27.16424 xmax: -76.68305 ymax: 35.66163
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(MAGR4.range) <- 4267 MAGR4_clipped = st_intersection(MAGR4.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = MAGR4_clipped, col = "red", linewidth = 2)+
geom_point(data = MAGR4.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/MAGR4.range.pdf", width = 12, height = 8)Subset for Magnolia macrophylla
MAMA2.occ = gbif %>%
filter(species == "Magnolia macrophylla") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
MAMA2.range = st_read("../../USTreeAtlas/shp/magnmacr/")## Reading layer `magnmacr' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/magnmacr'
## using driver `ESRI Shapefile'
## Simple feature collection with 25 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -92.80724 ymin: 30.45176 xmax: -78.51603 ymax: 39.02176
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(MAMA2.range) <- 4267 MAMA2_clipped = st_intersection(MAMA2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = MAMA2_clipped, col = "red", linewidth = 2)+
geom_point(data = MAMA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/MAMA2.range.pdf", width = 12, height = 8)Subset for Magnolia virginiana
MAVI2.occ = gbif %>%
filter(species == "Magnolia virginiana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
MAVI2.range = st_read("../../USTreeAtlas/shp/magnvirg/")## Reading layer `magnvirg' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/magnvirg'
## using driver `ESRI Shapefile'
## Simple feature collection with 82 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -95.6114 ymin: 25.71463 xmax: -72.75641 ymax: 40.91171
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(MAVI2.range) <- 4267 MAVI2_clipped = st_intersection(MAVI2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = MAVI2_clipped, col = "red", linewidth = 2)+
geom_point(data = MAVI2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/MAVI2.range.pdf", width = 12, height = 8)Subset for Morus rubra
MORU2.occ = gbif %>%
filter(species == "Morus rubra") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
MORU2.range = st_read("../../USTreeAtlas/shp/morurubr/")## Reading layer `morurubr' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/morurubr'
## using driver `ESRI Shapefile'
## Simple feature collection with 96 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -101.9412 ymin: 25.71206 xmax: -70.92392 ymax: 44.91562
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(MORU2.range) <- 4267 MORU2_clipped = st_intersection(MORU2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = MORU2_clipped, col = "red", linewidth = 2)+
geom_point(data = MORU2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/MORU2.range.pdf", width = 12, height = 8)Subset for Nyssa aquatica
NYAQ2.occ = gbif %>%
filter(species == "Nyssa aquatica") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
NYAQ2.range = st_read("../../USTreeAtlas/shp/nyssaqua/")## Reading layer `nyssaqua' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/nyssaqua'
## using driver `ESRI Shapefile'
## Simple feature collection with 48 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -95.18513 ymin: 29.33239 xmax: -75.44749 ymax: 39.09947
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(NYAQ2.range) <- 4267 NYAQ2_clipped = st_intersection(NYAQ2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = NYAQ2_clipped, col = "red", linewidth = 2)+
geom_point(data = NYAQ2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/NYAQ2.range.pdf", width = 12, height = 8)Subset for Nyssa biflora. No range map available.
NYBI.occ = gbif %>%
filter(species == "Nyssa biflora") %>%
select(species,decimalLatitude,decimalLongitude)ggplot()+
geom_sf(data = states.map, fill = "white")+
geom_point(data = NYBI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/NYBI.range.pdf", width = 12, height = 8)Subset for Nyssa ogeche
NYOG.occ = gbif %>%
filter(species == "Nyssa ogeche") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
NYOG.range = st_read("../../USTreeAtlas/shp/nyssogec/")## Reading layer `nyssogec' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/nyssogec'
## using driver `ESRI Shapefile'
## Simple feature collection with 8 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -86.51112 ymin: 29.52137 xmax: -80.47782 ymax: 32.7832
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(NYOG.range) <- 4267 NYOG_clipped = st_intersection(NYOG.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = NYOG_clipped, col = "red", linewidth = 2)+
geom_point(data = NYOG.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/NYOG.range.pdf", width = 12, height = 8)Subset for Nyssa sylvatica
NYSY.occ = gbif %>%
filter(species == "Nyssa sylvatica") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
NYSY.range = st_read("../../USTreeAtlas/shp/nysssylv/")## Reading layer `nysssylv' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/nysssylv'
## using driver `ESRI Shapefile'
## Simple feature collection with 98 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -98.92374 ymin: 16.64896 xmax: -69.58422 ymax: 44.80475
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(NYSY.range) <- 4267 NYSY_clipped = st_intersection(NYSY.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = NYSY_clipped, col = "red", linewidth = 2)+
geom_point(data = NYSY.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/NYSY.range.pdf", width = 12, height = 8)Subset for Ostrya virginiana
OSVI.occ = gbif %>%
filter(species == "Ostrya virginiana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
OSVI.range = st_read("../../USTreeAtlas/shp/ostrvirg/")## Reading layer `ostrvirg' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/ostrvirg'
## using driver `ESRI Shapefile'
## Simple feature collection with 145 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -109.4911 ymin: 13.70159 xmax: -59.68972 ymax: 50.35658
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(OSVI.range) <- 4267 OSVI_clipped = st_intersection(OSVI.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = OSVI_clipped, col = "red", linewidth = 2)+
geom_point(data = OSVI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/OSVI.range.pdf", width = 12, height = 8)Subset for Oxydendrum arboreum
OXAR.occ = gbif %>%
filter(species == "Oxydendrum arboreum") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
OXAR.range = st_read("../../USTreeAtlas/shp/oxydarbo/")## Reading layer `oxydarbo' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/oxydarbo'
## using driver `ESRI Shapefile'
## Simple feature collection with 34 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -94.92362 ymin: 29.90913 xmax: -75.31862 ymax: 40.52448
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(OXAR.range) <- 4267 OXAR_clipped = st_intersection(OXAR.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = OXAR_clipped, col = "red", linewidth = 2)+
geom_point(data = OXAR.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/OXAR.range.pdf", width = 12, height = 8)Subset for Paulownia tomentosa. No range map available.
PATO2.occ = gbif %>%
filter(species == "Paulownia tomentosa") %>%
select(species,decimalLatitude,decimalLongitude)ggplot()+
geom_sf(data = states.map, fill = "white")+
geom_point(data = PATO2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PATO2.range.pdf", width = 12, height = 8)Subset for Persea borbonia
PEBO.occ = gbif %>%
filter(species == "Persea borbonia") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PEBO.range = st_read("../../USTreeAtlas/shp/persborb/")## Reading layer `persborb' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/persborb'
## using driver `ESRI Shapefile'
## Simple feature collection with 94 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -97.96335 ymin: 24.6522 xmax: -75.1954 ymax: 38.55726
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PEBO.range) <- 4267 PEBO_clipped = st_intersection(PEBO.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PEBO_clipped, col = "red", linewidth = 2)+
geom_point(data = PEBO.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PEBO.range.pdf", width = 12, height = 8)Subset for Picea glauca
PIGL.occ = gbif %>%
filter(species == "Picea glauca") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PIGL.range = st_read("../../USTreeAtlas/shp/piceglau/")## Reading layer `piceglau' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/piceglau'
## using driver `ESRI Shapefile'
## Simple feature collection with 1449 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -163.9861 ymin: 42.71046 xmax: -52.63629 ymax: 69.57938
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PIGL.range) <- 4267 PIGL_clipped = st_intersection(PIGL.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PIGL_clipped, col = "red", linewidth = 2)+
geom_point(data = PIGL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PIGL.range.pdf", width = 12, height = 8)Subset for Picea mariana
PIMA.occ = gbif %>%
filter(species == "Picea mariana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PIMA.range = st_read("../../USTreeAtlas/shp/picemari/")## Reading layer `picemari' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/picemari'
## using driver `ESRI Shapefile'
## Simple feature collection with 1585 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -161.9129 ymin: 40.52763 xmax: -52.63629 ymax: 69.51231
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PIMA.range) <- 4267 PIMA_clipped = st_intersection(PIMA.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PIMA_clipped, col = "red", linewidth = 2)+
geom_point(data = PIMA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PIMA.range.pdf", width = 12, height = 8)Subset for Picea rubens
PIRU.occ = gbif %>%
filter(species == "Picea rubens") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PIRU.range = st_read("../../USTreeAtlas/shp/picerube/")## Reading layer `picerube' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/picerube'
## using driver `ESRI Shapefile'
## Simple feature collection with 326 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -83.7767 ymin: 35.23503 xmax: -59.69256 ymax: 48.01405
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PIRU.range) <- 4267 PIRU_clipped = st_intersection(PIRU.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PIRU_clipped, col = "red", linewidth = 2)+
geom_point(data = PIRU.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PIRU.range.pdf", width = 12, height = 8)Subset for Pinus banksiana
PIBA2.occ = gbif %>%
filter(species == "Pinus banksiana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PIBA2.range = st_read("../../USTreeAtlas/shp/pinubank/")## Reading layer `pinubank' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinubank'
## using driver `ESRI Shapefile'
## Simple feature collection with 395 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -127.2683 ymin: 41.47121 xmax: -60.3597 ymax: 65.29585
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PIBA2.range) <- 4267 PIBA2_clipped = st_intersection(PIBA2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PIBA2_clipped, col = "red", linewidth = 2)+
geom_point(data = PIBA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PIBA2.range.pdf", width = 12, height = 8)Subset for Pinus clausa
PICL.occ = gbif %>%
filter(species == "Pinus clausa") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PICL.range = st_read("../../USTreeAtlas/shp/pinuclau/")## Reading layer `pinuclau' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinuclau'
## using driver `ESRI Shapefile'
## Simple feature collection with 52 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -87.70842 ymin: 25.84471 xmax: -80.03255 ymax: 30.66881
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PICL.range) <- 4267 PICL_clipped = st_intersection(PICL.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PICL_clipped, col = "red", linewidth = 2)+
geom_point(data = PICL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PICL.range.pdf", width = 12, height = 8)Subset for Pinus echinata
PIEC2.occ = gbif %>%
filter(species == "Pinus echinata") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PIEC2.range = st_read("../../USTreeAtlas/shp/pinuechi/")## Reading layer `pinuechi' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinuechi'
## using driver `ESRI Shapefile'
## Simple feature collection with 68 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.62461 ymin: 29.79749 xmax: -73.2604 ymax: 41.11497
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PIEC2.range) <- 4267 PIEC2_clipped = st_intersection(PIEC2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PIEC2_clipped, col = "red", linewidth = 2)+
geom_point(data = PIEC2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PIEC2.range.pdf", width = 12, height = 8)Subset for Pinus elliottii
PIEL.occ = gbif %>%
filter(species == "Pinus elliottii") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PIEL.range = st_read("../../USTreeAtlas/shp/pinuelli/")## Reading layer `pinuelli' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinuelli'
## using driver `ESRI Shapefile'
## Simple feature collection with 129 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -90.77222 ymin: 24.58278 xmax: -79.60842 ymax: 33.28711
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PIEL.range) <- 4267 PIEL_clipped = st_intersection(PIEL.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PIEL_clipped, col = "red", linewidth = 2)+
geom_point(data = PIEL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PIEL.range.pdf", width = 12, height = 8)Subset for Pinus glabra
PIGL2.occ = gbif %>%
filter(species == "Pinus glabra") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PIGL2.range = st_read("../../USTreeAtlas/shp/pinuglab/")## Reading layer `pinuglab' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinuglab'
## using driver `ESRI Shapefile'
## Simple feature collection with 9 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -91.29132 ymin: 29.58212 xmax: -79.20943 ymax: 34.05908
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PIGL2.range) <- 4267 PIGL2_clipped = st_intersection(PIGL2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PIGL2_clipped, col = "red", linewidth = 2)+
geom_point(data = PIGL2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PIGL2.range.pdf", width = 12, height = 8)Subset for Pinus palustris
PIPA2.occ = gbif %>%
filter(species == "Pinus palustris") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PIPA2.range = st_read("../../USTreeAtlas/shp/pinupalu/")## Reading layer `pinupalu' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinupalu'
## using driver `ESRI Shapefile'
## Simple feature collection with 38 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -95.21806 ymin: 26.61921 xmax: -75.79999 ymax: 36.84856
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PIPA2.range) <- 4267 PIPA2_clipped = st_intersection(PIPA2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PIPA2_clipped, col = "red", linewidth = 2)+
geom_point(data = PIPA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PIPA2.range.pdf", width = 12, height = 8)Subset for Pinus pungens
PIPUS.occ = gbif %>%
filter(species == "Pinus pungens") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PIPUS.range = st_read("../../USTreeAtlas/shp/pinupung/")## Reading layer `pinupung' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinupung'
## using driver `ESRI Shapefile'
## Simple feature collection with 26 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -84.29883 ymin: 34.56113 xmax: -74.77565 ymax: 41.06224
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PIPUS.range) <- 4267 PIPUS_clipped = st_intersection(PIPUS.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PIPUS_clipped, col = "red", linewidth = 2)+
geom_point(data = PIPUS.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PIPUS.range.pdf", width = 12, height = 8)Subset for Pinus resinosa
PIRE.occ = gbif %>%
filter(species == "Pinus resinosa") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PIRE.range = st_read("../../USTreeAtlas/shp/pinuresi/")## Reading layer `pinuresi' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinuresi'
## using driver `ESRI Shapefile'
## Simple feature collection with 485 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.63834 ymin: 38.64299 xmax: -54.26736 ymax: 51.25153
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PIRE.range) <- 4267 PIRE_clipped = st_intersection(PIRE.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PIRE_clipped, col = "red", linewidth = 2)+
geom_point(data = PIRE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PIRE.range.pdf", width = 12, height = 8)Subset for Pinus rigida
PIRI.occ = gbif %>%
filter(species == "Pinus rigida") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PIRI.range = st_read("../../USTreeAtlas/shp/pinurigi/")## Reading layer `pinurigi' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinurigi'
## using driver `ESRI Shapefile'
## Simple feature collection with 154 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -88.39619 ymin: 33.68433 xmax: -68.12421 ymax: 45.25713
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PIRI.range) <- 4267 PIRI_clipped = st_intersection(PIRI.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PIRI_clipped, col = "red", linewidth = 2)+
geom_point(data = PIRI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PIRI.range.pdf", width = 12, height = 8)Subset for Pinus serotina
PISE.occ = gbif %>%
filter(species == "Pinus serotina") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PISE.range = st_read("../../USTreeAtlas/shp/pinusero/")## Reading layer `pinusero' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinusero'
## using driver `ESRI Shapefile'
## Simple feature collection with 40 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -87.52861 ymin: 28.02775 xmax: -74.69576 ymax: 39.87193
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PISE.range) <- 4267 PISE_clipped = st_intersection(PISE.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PISE_clipped, col = "red", linewidth = 2)+
geom_point(data = PISE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PISE.range.pdf", width = 12, height = 8)Subset for Pinus strobus
PIST.occ = gbif %>%
filter(species == "Pinus strobus") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PIST.range = st_read("../../USTreeAtlas/shp/pinustrb/")## Reading layer `pinustrb' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinustrb'
## using driver `ESRI Shapefile'
## Simple feature collection with 264 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -95.83807 ymin: 34.27946 xmax: -52.61445 ymax: 51.11618
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PIST.range) <- 4267 PIST_clipped = st_intersection(PIST.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PIST_clipped, col = "red", linewidth = 2)+
geom_point(data = PIST.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PIST.range.pdf", width = 12, height = 8)Subset for Pinus taeda
PITA.occ = gbif %>%
filter(species == "Pinus taeda") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PITA.range = st_read("../../USTreeAtlas/shp/pinutaed/")## Reading layer `pinutaed' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinutaed'
## using driver `ESRI Shapefile'
## Simple feature collection with 58 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -97.6138 ymin: 28.22558 xmax: -75.01347 ymax: 39.50686
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PITA.range) <- 4267 PITA_clipped = st_intersection(PITA.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PITA_clipped, col = "red", linewidth = 2)+
geom_point(data = PITA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PITA.range.pdf", width = 12, height = 8)Subset for Pinus virginiana
PIVI2.occ = gbif %>%
filter(species == "Pinus virginiana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PIVI2.range = st_read("../../USTreeAtlas/shp/pinuvirg/")## Reading layer `pinuvirg' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinuvirg'
## using driver `ESRI Shapefile'
## Simple feature collection with 36 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -89.41452 ymin: 32.30019 xmax: -73.15847 ymax: 41.37672
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PIVI2.range) <- 4267 PIVI2_clipped = st_intersection(PIVI2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PIVI2_clipped, col = "red", linewidth = 2)+
geom_point(data = PIVI2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PIVI2.range.pdf", width = 12, height = 8)Subset for Planera aquatica
PLAQ.occ = gbif %>%
filter(species == "Planera aquatica") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PLAQ.range = st_read("../../USTreeAtlas/shp/planaqua/")## Reading layer `planaqua' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/planaqua'
## using driver `ESRI Shapefile'
## Simple feature collection with 17 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.84963 ymin: 29.49777 xmax: -78.27449 ymax: 37.5579
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PLAQ.range) <- 4267 PLAQ_clipped = st_intersection(PLAQ.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PLAQ_clipped, col = "red", linewidth = 2)+
geom_point(data = PLAQ.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PLAQ.range.pdf", width = 12, height = 8)Subset for Platanus occidentalis
PLOC.occ = gbif %>%
filter(species == "Platanus occidentalis") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PLOC.range = st_read("../../USTreeAtlas/shp/platocci/")## Reading layer `platocci' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/platocci'
## using driver `ESRI Shapefile'
## Simple feature collection with 77 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -102.2446 ymin: 21.80153 xmax: -69.63906 ymax: 44.61452
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PLOC.range) <- 4267 PLOC_clipped = st_intersection(PLOC.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PLOC_clipped, col = "red", linewidth = 2)+
geom_point(data = PLOC.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PLOC.range.pdf", width = 12, height = 8)Subset for Populus balsamifera
POBA2.occ = gbif %>%
filter(species == "Populus balsamifera") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
POBA2.range = st_read("../../USTreeAtlas/shp/popubals/")## Reading layer `popubals' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/popubals'
## using driver `ESRI Shapefile'
## Simple feature collection with 415 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -163.9213 ymin: 38.58866 xmax: -52.61445 ymax: 69.24193
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(POBA2.range) <- 4267 POBA2_clipped = st_intersection(POBA2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = POBA2_clipped, col = "red", linewidth = 2)+
geom_point(data = POBA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/POBA2.range.pdf", width = 12, height = 8)Subset for Populus deltoides
PODE3.occ = gbif %>%
filter(species == "Populus deltoides") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PODE3.range = st_read("../../USTreeAtlas/shp/popudelt/")## Reading layer `popudelt' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/popudelt'
## using driver `ESRI Shapefile'
## Simple feature collection with 79 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -114.3961 ymin: 27.32564 xmax: -72.24688 ymax: 52.14172
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PODE3.range) <- 4267 PODE3_clipped = st_intersection(PODE3.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PODE3_clipped, col = "red", linewidth = 2)+
geom_point(data = PODE3.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PODE3.range.pdf", width = 12, height = 8)Subset for Populus grandidentata
POGR4.occ = gbif %>%
filter(species == "Populus grandidentata") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
POGR4.range = st_read("../../USTreeAtlas/shp/popugran/")## Reading layer `popugran' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/popugran'
## using driver `ESRI Shapefile'
## Simple feature collection with 157 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.50261 ymin: 35.13372 xmax: -59.68972 ymax: 49.99322
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(POGR4.range) <- 4267 POGR4_clipped = st_intersection(POGR4.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = POGR4_clipped, col = "red", linewidth = 2)+
geom_point(data = POGR4.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/POGR4.range.pdf", width = 12, height = 8)Subset for Populus heterophylla
POHE4.occ = gbif %>%
filter(species == "Populus heterophylla") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
POHE4.range = st_read("../../USTreeAtlas/shp/popuhete/")## Reading layer `popuhete' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/popuhete'
## using driver `ESRI Shapefile'
## Simple feature collection with 79 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -92.28237 ymin: 29.70388 xmax: -71.85611 ymax: 42.36729
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(POHE4.range) <- 4267 POHE4_clipped = st_intersection(POHE4.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = POHE4_clipped, col = "red", linewidth = 2)+
geom_point(data = POHE4.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/POHE4.range.pdf", width = 12, height = 8)Subset for Populus tremuloides
POTR5.occ = gbif %>%
filter(species == "Populus tremuloides") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
POTR5.range = st_read("../../USTreeAtlas/shp/poputrem/")## Reading layer `poputrem' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/poputrem'
## using driver `ESRI Shapefile'
## Simple feature collection with 534 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -163.7564 ymin: 20.77411 xmax: -52.61445 ymax: 69.2738
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(POTR5.range) <- 4267 POTR5_clipped = st_intersection(POTR5.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = POTR5_clipped, col = "red", linewidth = 2)+
geom_point(data = POTR5.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/POTR5.range.pdf", width = 12, height = 8)Subset for Prunus americana
PRAM.occ = gbif %>%
filter(species == "Prunus americana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PRAM.range = st_read("../../USTreeAtlas/shp/prunamer/")## Reading layer `prunamer' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/prunamer'
## using driver `ESRI Shapefile'
## Simple feature collection with 71 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -114.464 ymin: 28.57408 xmax: -70.17076 ymax: 50.03462
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PRAM.range) <- 4267 PRAM_clipped = st_intersection(PRAM.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PRAM_clipped, col = "red", linewidth = 2)+
geom_point(data = PRAM.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PRAM.range.pdf", width = 12, height = 8)Subset for Prunus pensylvanica
PRPE2.occ = gbif %>%
filter(species == "Prunus pensylvanica") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PRPE2.range = st_read("../../USTreeAtlas/shp/prunpens/")## Reading layer `prunpens' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/prunpens'
## using driver `ESRI Shapefile'
## Simple feature collection with 336 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -128.2946 ymin: 34.80802 xmax: -52.61445 ymax: 62.31179
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PRPE2.range) <- 4267 PRPE2_clipped = st_intersection(PRPE2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PRPE2_clipped, col = "red", linewidth = 2)+
geom_point(data = PRPE2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PRPE2.range.pdf", width = 12, height = 8)Subset for Prunus serotina
PRSE2.occ = gbif %>%
filter(species == "Prunus serotina") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PRSE2.range = st_read("../../USTreeAtlas/shp/prunsero/")## Reading layer `prunsero' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/prunsero'
## using driver `ESRI Shapefile'
## Simple feature collection with 200 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -113.9141 ymin: 14.56384 xmax: -63.3591 ymax: 47.8594
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PRSE2.range) <- 4267 PRSE2_clipped = st_intersection(PRSE2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PRSE2_clipped, col = "red", linewidth = 2)+
geom_point(data = PRSE2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PRSE2.range.pdf", width = 12, height = 8)Subset for Prunus virginiana
PRVI.occ = gbif %>%
filter(species == "Prunus virginiana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
PRVI.range = st_read("../../USTreeAtlas/shp/prunvirg/")## Reading layer `prunvirg' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/prunvirg'
## using driver `ESRI Shapefile'
## Simple feature collection with 397 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -130.5808 ymin: 30.24095 xmax: -53.03568 ymax: 61.54309
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(PRVI.range) <- 4267 PRVI_clipped = st_intersection(PRVI.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = PRVI_clipped, col = "red", linewidth = 2)+
geom_point(data = PRVI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/PRVI.range.pdf", width = 12, height = 8)Subset for Quercus alba
QUAL.occ = gbif %>%
filter(species == "Quercus alba") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUAL.range = st_read("../../USTreeAtlas/shp/queralba/")## Reading layer `queralba' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/queralba'
## using driver `ESRI Shapefile'
## Simple feature collection with 82 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.34524 ymin: 29.63241 xmax: -69.11673 ymax: 46.45134
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUAL.range) <- 4267 QUAL_clipped = st_intersection(QUAL.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUAL_clipped, col = "red", linewidth = 2)+
geom_point(data = QUAL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUAL.range.pdf", width = 12, height = 8)Subset for Quercus bicolor
QUBI.occ = gbif %>%
filter(species == "Quercus bicolor") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUBI.range = st_read("../../USTreeAtlas/shp/querbico/")## Reading layer `querbico' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querbico'
## using driver `ESRI Shapefile'
## Simple feature collection with 49 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.39155 ymin: 35.17383 xmax: -70.02106 ymax: 46.36058
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUBI.range) <- 4267 QUBI_clipped = st_intersection(QUBI.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUBI_clipped, col = "red", linewidth = 2)+
geom_point(data = QUBI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUBI.range.pdf", width = 12, height = 8)Subset for Quercus coccinea
QUCO2.occ = gbif %>%
filter(species == "Quercus coccinea") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUCO2.range = st_read("../../USTreeAtlas/shp/quercocc/")## Reading layer `quercocc' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quercocc'
## using driver `ESRI Shapefile'
## Simple feature collection with 76 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -92.35469 ymin: 31.22379 xmax: -69.87 ymax: 44.1226
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUCO2.range) <- 4267 QUCO2_clipped = st_intersection(QUCO2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUCO2_clipped, col = "red", linewidth = 2)+
geom_point(data = QUCO2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUCO2.range.pdf", width = 12, height = 8)Subset for Quercus ellipsoidalis
QUEL.occ = gbif %>%
filter(species == "Quercus ellipsoidalis") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUEL.range = st_read("../../USTreeAtlas/shp/querelli/")## Reading layer `querelli' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querelli'
## using driver `ESRI Shapefile'
## Simple feature collection with 22 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.91765 ymin: 40.28243 xmax: -82.88901 ymax: 48.66541
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUEL.range) <- 4267 QUEL_clipped = st_intersection(QUEL.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUEL_clipped, col = "red", linewidth = 2)+
geom_point(data = QUEL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUEL.range.pdf", width = 12, height = 8)Subset for Quercus falcata
QUFA.occ = gbif %>%
filter(species == "Quercus falcata") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUFA.range = st_read("../../USTreeAtlas/shp/querfalc/")## Reading layer `querfalc' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querfalc'
## using driver `ESRI Shapefile'
## Simple feature collection with 58 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.82243 ymin: 28.8895 xmax: -73.46515 ymax: 40.73988
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUFA.range) <- 4267 QUFA_clipped = st_intersection(QUFA.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUFA_clipped, col = "red", linewidth = 2)+
geom_point(data = QUFA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUFA.range.pdf", width = 12, height = 8)Subset for Quercus ilicifolia
QUIL.occ = gbif %>%
filter(species == "Quercus ilicifolia") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUIL.range = st_read("../../USTreeAtlas/shp/querilic/")## Reading layer `querilic' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querilic'
## using driver `ESRI Shapefile'
## Simple feature collection with 46 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -81.30976 ymin: 35.22381 xmax: -67.91949 ymax: 44.84274
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUIL.range) <- 4267 QUIL_clipped = st_intersection(QUIL.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUIL_clipped, col = "red", linewidth = 2)+
geom_point(data = QUIL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUIL.range.pdf", width = 12, height = 8)Subset for Quercus imbricaria
QUIM.occ = gbif %>%
filter(species == "Quercus imbricaria") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUIM.range = st_read("../../USTreeAtlas/shp/querimbr/")## Reading layer `querimbr' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querimbr'
## using driver `ESRI Shapefile'
## Simple feature collection with 38 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.96888 ymin: 31.03043 xmax: -75.56732 ymax: 42.31751
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUIM.range) <- 4267 QUIM_clipped = st_intersection(QUIM.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUIM_clipped, col = "red", linewidth = 2)+
geom_point(data = QUIM.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUIM.range.pdf", width = 12, height = 8)Subset for Quercus incana
QUIN.occ = gbif %>%
filter(species == "Quercus incana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUIN.range = st_read("../../USTreeAtlas/shp/querinca/")## Reading layer `querinca' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querinca'
## using driver `ESRI Shapefile'
## Simple feature collection with 55 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -97.78163 ymin: 26.41639 xmax: -75.77622 ymax: 36.8999
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUIN.range) <- 4267 QUIN_clipped = st_intersection(QUIN.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUIN_clipped, col = "red", linewidth = 2)+
geom_point(data = QUIN.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUIN.range.pdf", width = 12, height = 8)Subset for Quercus laevis
QULA2.occ = gbif %>%
filter(species == "Quercus laevis") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QULA2.range = st_read("../../USTreeAtlas/shp/querlaev/")## Reading layer `querlaev' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querlaev'
## using driver `ESRI Shapefile'
## Simple feature collection with 44 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -90.19388 ymin: 26.39276 xmax: -75.75282 ymax: 36.97728
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QULA2.range) <- 4267 QULA2_clipped = st_intersection(QULA2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QULA2_clipped, col = "red", linewidth = 2)+
geom_point(data = QULA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QULA2.range.pdf", width = 12, height = 8)Subset for Quercus laurifolia
QULA3.occ = gbif %>%
filter(species == "Quercus laurifolia") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QULA3.range = st_read("../../USTreeAtlas/shp/querlaur/")## Reading layer `querlaur' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querlaur'
## using driver `ESRI Shapefile'
## Simple feature collection with 73 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -97.66097 ymin: 25.38946 xmax: -75.44749 ymax: 37.41531
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QULA3.range) <- 4267 QULA3_clipped = st_intersection(QULA3.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QULA3_clipped, col = "red", linewidth = 2)+
geom_point(data = QULA3.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QULA3.range.pdf", width = 12, height = 8)Subset for Quercus lyrata
QULY.occ = gbif %>%
filter(species == "Quercus lyrata") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QULY.range = st_read("../../USTreeAtlas/shp/querlyra/")## Reading layer `querlyra' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querlyra'
## using driver `ESRI Shapefile'
## Simple feature collection with 41 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -97.41548 ymin: 28.61861 xmax: -75.03833 ymax: 41.84388
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QULY.range) <- 4267 QULY_clipped = st_intersection(QULY.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QULY_clipped, col = "red", linewidth = 2)+
geom_point(data = QULY.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QULY.range.pdf", width = 12, height = 8)Subset for Quercus macrocarpa
QUMA2.occ = gbif %>%
filter(species == "Quercus macrocarpa") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUMA2.range = st_read("../../USTreeAtlas/shp/quermacr/")## Reading layer `quermacr' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quermacr'
## using driver `ESRI Shapefile'
## Simple feature collection with 41 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -104.4497 ymin: 28.01439 xmax: -66.08747 ymax: 52.74889
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUMA2.range) <- 4267 QUMA2_clipped = st_intersection(QUMA2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUMA2_clipped, col = "red", linewidth = 2)+
geom_point(data = QUMA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUMA2.range.pdf", width = 12, height = 8)Subset for Quercus margaretta. No range map available
QUMA6.occ = gbif %>%
filter(species == "Quercus margaretta") %>%
select(species,decimalLatitude,decimalLongitude)ggplot()+
geom_sf(data = states.map, fill = "white")+
geom_point(data = QUMA6.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUMA6.range.pdf", width = 12, height = 8)Subset for Quercus marilandica
QUMA3.occ = gbif %>%
filter(species == "Quercus marilandica") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUMA3.range = st_read("../../USTreeAtlas/shp/quermari/")## Reading layer `quermari' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quermari'
## using driver `ESRI Shapefile'
## Simple feature collection with 53 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -99.95763 ymin: 28.12027 xmax: -73.53086 ymax: 42.7558
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUMA3.range) <- 4267 QUMA3_clipped = st_intersection(QUMA3.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUMA3_clipped, col = "red", linewidth = 2)+
geom_point(data = QUMA3.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUMA3.range.pdf", width = 12, height = 8)Subset for Quercus michauxii (gbif) and Quercus prinus
QUMI.occ = gbif %>%
filter(species == "Quercus michauxii") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUMI.range = st_read("../../USTreeAtlas/shp/quermich/")## Reading layer `quermich' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quermich'
## using driver `ESRI Shapefile'
## Simple feature collection with 64 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -95.4984 ymin: 28.78528 xmax: -74.28231 ymax: 40.95467
## CRS: NAQUPR2.range = st_read("../../USTreeAtlas/shp/querprin/")## Reading layer `querprin' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querprin'
## using driver `ESRI Shapefile'
## Simple feature collection with 29 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -89.98369 ymin: 31.95997 xmax: -70.50323 ymax: 44.57826
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUMI.range) <- 4267
st_crs(QUPR2.range) <- 4267 QUMI_clipped = st_intersection(QUMI.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesQUPR2_clipped = st_intersection(QUPR2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUMI_clipped, col = "red", linewidth = 2)+
geom_sf(data = QUPR2_clipped, col = "red", linewidth = 2)+
geom_point(data = QUMI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUMI.range.pdf", width = 12, height = 8)Subset for Quercus minima. No range map available.
QUMI2.occ = gbif %>%
filter(species == "Quercus minima") %>%
select(species,decimalLatitude,decimalLongitude)ggplot()+
geom_sf(data = states.map, fill = "white")+
geom_point(data = QUMI2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUMI2.range.pdf", width = 12, height = 8)Subset for Quercus muehlenbergii
QUMU.occ = gbif %>%
filter(species == "Quercus muehlenbergii") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUMU.range = st_read("../../USTreeAtlas/shp/quermueh/")## Reading layer `quermueh' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quermueh'
## using driver `ESRI Shapefile'
## Simple feature collection with 32 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -105.2417 ymin: 24.78449 xmax: -72.23889 ymax: 44.69553
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUMU.range) <- 4267 QUMU_clipped = st_intersection(QUMU.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUMU_clipped, col = "red", linewidth = 2)+
geom_point(data = QUMU.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUMU.range.pdf", width = 12, height = 8)Subset for Quercus nigra
QUNI.occ = gbif %>%
filter(species == "Quercus nigra") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUNI.range = st_read("../../USTreeAtlas/shp/quernigr/")## Reading layer `quernigr' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quernigr'
## using driver `ESRI Shapefile'
## Simple feature collection with 45 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -97.70858 ymin: 27.65694 xmax: -74.65406 ymax: 39.40351
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUNI.range) <- 4267 QUNI_clipped = st_intersection(QUNI.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUNI_clipped, col = "red", linewidth = 2)+
geom_point(data = QUNI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUNI.range.pdf", width = 12, height = 8)Subset for Quercus pagoda. No range map available.
QUPA5.occ = gbif %>%
filter(species == "Quercus pagoda") %>%
select(species,decimalLatitude,decimalLongitude)ggplot()+
geom_sf(data = states.map, fill = "white")+
geom_point(data = QUPA5.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUPA5.range.pdf", width = 12, height = 8)Subset for Quercus palustris
QUPA2.occ = gbif %>%
filter(species == "Quercus palustris") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUPA2.range = st_read("../../USTreeAtlas/shp/querpalu/")## Reading layer `querpalu' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querpalu'
## using driver `ESRI Shapefile'
## Simple feature collection with 31 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.838 ymin: 34.436 xmax: -71.4752 ymax: 43.71538
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUPA2.range) <- 4267 QUPA2_clipped = st_intersection(QUPA2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUPA2_clipped, col = "red", linewidth = 2)+
geom_point(data = QUPA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUPA2.range.pdf", width = 12, height = 8)Subset for Quercus phellos
QUPH.occ = gbif %>%
filter(species == "Quercus phellos") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUPH.range = st_read("../../USTreeAtlas/shp/querphel/")## Reading layer `querphel' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querphel'
## using driver `ESRI Shapefile'
## Simple feature collection with 42 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.26867 ymin: 28.91285 xmax: -73.14029 ymax: 40.74458
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUPH.range) <- 4267 QUPH_clipped = st_intersection(QUPH.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUPH_clipped, col = "red", linewidth = 2)+
geom_point(data = QUPH.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUPH.range.pdf", width = 12, height = 8)Subset for Quercus rubra
QURU.occ = gbif %>%
filter(species == "Quercus rubra") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QURU.range = st_read("../../USTreeAtlas/shp/querrubr/")## Reading layer `querrubr' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querrubr'
## using driver `ESRI Shapefile'
## Simple feature collection with 155 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -97.78265 ymin: 30.7561 xmax: -59.68972 ymax: 49.08564
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QURU.range) <- 4267 QURU_clipped = st_intersection(QURU.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QURU_clipped, col = "red", linewidth = 2)+
geom_point(data = QURU.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QURU.range.pdf", width = 12, height = 8)Subset for Quercus shumardii
QUSH.occ = gbif %>%
filter(species == "Quercus shumardii") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUSH.range = st_read("../../USTreeAtlas/shp/quershum/")## Reading layer `quershum' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quershum'
## using driver `ESRI Shapefile'
## Simple feature collection with 47 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -100.2109 ymin: 28.89789 xmax: -76.45997 ymax: 42.35025
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUSH.range) <- 4267 QUSH_clipped = st_intersection(QUSH.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUSH_clipped, col = "red", linewidth = 2)+
geom_point(data = QUSH.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUSH.range.pdf", width = 12, height = 8)Subset for Quercus sinuata. No range map available.
QUSIS.occ = gbif %>%
filter(species == "Quercus sinuata") %>%
select(species,decimalLatitude,decimalLongitude)ggplot()+
geom_sf(data = states.map, fill = "white")+
geom_point(data = QUSIS.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUSIS.range.pdf", width = 12, height = 8)Subset for Quercus stellata
QUST.occ = gbif %>%
filter(species == "Quercus stellata") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUST.range = st_read("../../USTreeAtlas/shp/querstel/")## Reading layer `querstel' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querstel'
## using driver `ESRI Shapefile'
## Simple feature collection with 74 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -101.3744 ymin: 27.61028 xmax: -69.95983 ymax: 41.76135
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUST.range) <- 4267 QUST_clipped = st_intersection(QUST.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUST_clipped, col = "red", linewidth = 2)+
geom_point(data = QUST.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUST.range.pdf", width = 12, height = 8)Subset for Quercus texana. No range map available.
QUTE.occ = gbif %>%
filter(species == "Quercus texana") %>%
select(species,decimalLatitude,decimalLongitude)ggplot()+
geom_sf(data = states.map, fill = "white")+
geom_point(data = QUTE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUTE.range.pdf", width = 12, height = 8)Subset for Quercus velutina
QUVE.occ = gbif %>%
filter(species == "Quercus velutina") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUVE.range = st_read("../../USTreeAtlas/shp/quervelu/")## Reading layer `quervelu' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quervelu'
## using driver `ESRI Shapefile'
## Simple feature collection with 72 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -97.28247 ymin: 29.83254 xmax: -69.37499 ymax: 45.49608
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUVE.range) <- 4267 QUVE_clipped = st_intersection(QUVE.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUVE_clipped, col = "red", linewidth = 2)+
geom_point(data = QUVE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUVE.range.pdf", width = 12, height = 8)Subset for Quercus virginiana
QUVI.occ = gbif %>%
filter(species == "Quercus virginiana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
QUVI.range = st_read("../../USTreeAtlas/shp/quervirg/")## Reading layer `quervirg' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quervirg'
## using driver `ESRI Shapefile'
## Simple feature collection with 94 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -102.121 ymin: 22.41014 xmax: -75.71973 ymax: 37.46014
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(QUVI.range) <- 4267 QUVI_clipped = st_intersection(QUVI.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = QUVI_clipped, col = "red", linewidth = 2)+
geom_point(data = QUVI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/QUVI.range.pdf", width = 12, height = 8)Subset for Robinia pseudoacacia
ROPS.occ = gbif %>%
filter(species == "Robinia pseudoacacia") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ROPS.range = st_read("../../USTreeAtlas/shp/robipseu/")## Reading layer `robipseu' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/robipseu'
## using driver `ESRI Shapefile'
## Simple feature collection with 12 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -95.69039 ymin: 32.40928 xmax: -76.87229 ymax: 41.12437
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ROPS.range) <- 4267 ROPS_clipped = st_intersection(ROPS.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ROPS_clipped, col = "red", linewidth = 2)+
geom_point(data = ROPS.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ROPS.range.pdf", width = 12, height = 8)Subset for Salix caroliniana
CACA5.occ = gbif %>%
filter(species == "Salix caroliniana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
CACA5.range = st_read("../../USTreeAtlas/shp/salicaro/")## Reading layer `salicaro' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/salicaro'
## using driver `ESRI Shapefile'
## Simple feature collection with 100 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -100.292 ymin: 25.10583 xmax: -75.52333 ymax: 40.99319
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(CACA5.range) <- 4267 CACA5_clipped = st_intersection(CACA5.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = CACA5_clipped, col = "red", linewidth = 2)+
geom_point(data = CACA5.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/CACA5.range.pdf", width = 12, height = 8)Subset for Salix nigra
SANI.occ = gbif %>%
filter(species == "Salix nigra") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
SANI.range = st_read("../../USTreeAtlas/shp/salinigr/")## Reading layer `salinigr' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/salinigr'
## using driver `ESRI Shapefile'
## Simple feature collection with 151 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -122.7956 ymin: 23.5848 xmax: -65.17285 ymax: 46.61042
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(SANI.range) <- 4267 SANI_clipped = st_intersection(SANI.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = SANI_clipped, col = "red", linewidth = 2)+
geom_point(data = SANI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/SANI.range.pdf", width = 12, height = 8)Subset for Sassafras albidum
SAAL5.occ = gbif %>%
filter(species == "Sassafras albidum") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
SAAL5.range = st_read("../../USTreeAtlas/shp/sassalbi/")## Reading layer `sassalbi' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/sassalbi'
## using driver `ESRI Shapefile'
## Simple feature collection with 74 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -96.55696 ymin: 28.49213 xmax: -69.92722 ymax: 44.81874
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(SAAL5.range) <- 4267 SAAL5_clipped = st_intersection(SAAL5.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = SAAL5_clipped, col = "red", linewidth = 2)+
geom_point(data = SAAL5.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/SAAL5.range.pdf", width = 12, height = 8)Subset for Sorbus americana
SOAM3.occ = gbif %>%
filter(species == "Sorbus americana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
SOAM3.range = st_read("../../USTreeAtlas/shp/sorbamer/")## Reading layer `sorbamer' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/sorbamer'
## using driver `ESRI Shapefile'
## Simple feature collection with 288 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -94.15333 ymin: 34.78817 xmax: -52.61445 ymax: 52.65512
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(SOAM3.range) <- 4267 SOAM3_clipped = st_intersection(SOAM3.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = SOAM3_clipped, col = "red", linewidth = 2)+
geom_point(data = SOAM3.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/SOAM3.range.pdf", width = 12, height = 8)Subset for Taxodium distichum (gbif) and Taxodium ascendens
TADI2.occ = gbif %>%
filter(species == "Taxodium distichum") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
TADI2.range = st_read("../../USTreeAtlas/shp/taxodist/")## Reading layer `taxodist' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/taxodist'
## using driver `ESRI Shapefile'
## Simple feature collection with 85 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -99.82847 ymin: 24.54479 xmax: -75.04254 ymax: 38.92478
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(TADI2.range) <- 4267 TADI2_clipped = st_intersection(TADI2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = TADI2_clipped, col = "red", linewidth = 2)+
geom_point(data = TADI2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/TADI2.range.pdf", width = 12, height = 8)Subset for Thuja occidentalis
THOC2.occ = gbif %>%
filter(species == "Thuja occidentalis") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
THOC2.range = st_read("../../USTreeAtlas/shp/thujocci/")## Reading layer `thujocci' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/thujocci'
## using driver `ESRI Shapefile'
## Simple feature collection with 134 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -102.1047 ymin: 35.59812 xmax: -63.26956 ymax: 55.0168
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(THOC2.range) <- 4267 THOC2_clipped = st_intersection(THOC2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = THOC2_clipped, col = "red", linewidth = 2)+
geom_point(data = THOC2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/THOC2.range.pdf", width = 12, height = 8)Subset for Tilia americana
TIAM.occ = gbif %>%
filter(species == "Tilia americana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
TIAM.range = st_read("../../USTreeAtlas/shp/tiliamer/")## Reading layer `tiliamer' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/tiliamer'
## using driver `ESRI Shapefile'
## Simple feature collection with 94 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -100.9682 ymin: 34.97041 xmax: -65.77465 ymax: 49.99133
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(TIAM.range) <- 4267 TIAM_clipped = st_intersection(TIAM.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = TIAM_clipped, col = "red", linewidth = 2)+
geom_point(data = TIAM.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/TIAM.range.pdf", width = 12, height = 8)Subset for Triadica sebifera. No range map available.
TRSE6.occ = gbif %>%
filter(species == "Triadica sebifera") %>%
select(species,decimalLatitude,decimalLongitude)ggplot()+
geom_sf(data = states.map, fill = "white")+
geom_point(data = TRSE6.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/TRSE6.range.pdf", width = 12, height = 8)Subset for Tsuga canadensis
TSCA.occ = gbif %>%
filter(species == "Tsuga canadensis") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
TSCA.range = st_read("../../USTreeAtlas/shp/tsugcana/")## Reading layer `tsugcana' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/tsugcana'
## using driver `ESRI Shapefile'
## Simple feature collection with 180 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -93.74585 ymin: 33.5981 xmax: -59.68972 ymax: 48.33391
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(TSCA.range) <- 4267 TSCA_clipped = st_intersection(TSCA.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = TSCA_clipped, col = "red", linewidth = 2)+
geom_point(data = TSCA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/TSCA.range.pdf", width = 12, height = 8)Subset for Tsuga caroliniana
TSCA2.occ = gbif %>%
filter(species == "Tsuga caroliniana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
TSCA2.range = st_read("../../USTreeAtlas/shp/tsugcaro/")## Reading layer `tsugcaro' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/tsugcaro'
## using driver `ESRI Shapefile'
## Simple feature collection with 9 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -83.89653 ymin: 34.66754 xmax: -79.44339 ymax: 37.66805
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(TSCA2.range) <- 4267 TSCA2_clipped = st_intersection(TSCA2.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = TSCA2_clipped, col = "red", linewidth = 2)+
geom_point(data = TSCA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/TSCA2.range.pdf", width = 12, height = 8)Subset for Ulmus alata
ULAL.occ = gbif %>%
filter(species == "Ulmus alata") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ULAL.range = st_read("../../USTreeAtlas/shp/ulmualat/")## Reading layer `ulmualat' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/ulmualat'
## using driver `ESRI Shapefile'
## Simple feature collection with 38 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -98.37525 ymin: 27.26967 xmax: -75.44749 ymax: 39.47826
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ULAL.range) <- 4267 ULAL_clipped = st_intersection(ULAL.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ULAL_clipped, col = "red", linewidth = 2)+
geom_point(data = ULAL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ULAL.range.pdf", width = 12, height = 8)Subset for Ulmus americana
ULAM.occ = gbif %>%
filter(species == "Ulmus americana") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ULAM.range = st_read("../../USTreeAtlas/shp/ulmuamer/")## Reading layer `ulmuamer' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/ulmuamer'
## using driver `ESRI Shapefile'
## Simple feature collection with 207 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -106.7759 ymin: 26.73745 xmax: -59.68972 ymax: 53.81448
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ULAM.range) <- 4267 ULAM_clipped = st_intersection(ULAM.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ULAM_clipped, col = "red", linewidth = 2)+
geom_point(data = ULAM.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ULAM.range.pdf", width = 12, height = 8)Subset for Ulmus crassifolia
ULCR.occ = gbif %>%
filter(species == "Ulmus crassifolia") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ULCR.range = st_read("../../USTreeAtlas/shp/ulmucras/")## Reading layer `ulmucras' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/ulmucras'
## using driver `ESRI Shapefile'
## Simple feature collection with 15 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -100.5271 ymin: 25.47785 xmax: -89.93496 ymax: 36.42969
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ULCR.range) <- 4267 ULCR_clipped = st_intersection(ULCR.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ULCR_clipped, col = "red", linewidth = 2)+
geom_point(data = ULCR.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ULCR.range.pdf", width = 12, height = 8)Subset for Ulmus rubra
ULRU.occ = gbif %>%
filter(species == "Ulmus rubra") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ULRU.range = st_read("../../USTreeAtlas/shp/ulmurubr/")## Reading layer `ulmurubr' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/ulmurubr'
## using driver `ESRI Shapefile'
## Simple feature collection with 50 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -101.2964 ymin: 29.64866 xmax: -68.63799 ymax: 48.01849
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ULRU.range) <- 4267 ULRU_clipped = st_intersection(ULRU.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ULRU_clipped, col = "red", linewidth = 2)+
geom_point(data = ULRU.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ULRU.range.pdf", width = 12, height = 8)Subset for Ulmus thomasii
ULTH.occ = gbif %>%
filter(species == "Ulmus thomasii") %>%
select(species,decimalLatitude,decimalLongitude)Read in range map
ULTH.range = st_read("../../USTreeAtlas/shp/ulmuthom/")## Reading layer `ulmuthom' from data source
## `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/ulmuthom'
## using driver `ESRI Shapefile'
## Simple feature collection with 49 features and 5 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: -99.79724 ymin: 35.06176 xmax: -72.06438 ymax: 47.69813
## CRS: NAThe projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters
st_crs(ULTH.range) <- 4267 ULTH_clipped = st_intersection(ULTH.range, states.map)## Warning: attribute variables are assumed to be spatially constant throughout
## all geometriesggplot()+
geom_sf(data = states.map, fill = "white")+
geom_sf(data = ULTH_clipped, col = "red", linewidth = 2)+
geom_point(data = ULTH.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
theme_classic()
ggsave("../Plots/ULTH.range.pdf", width = 12, height = 8)